Siegert pseudostate perturbation theory: one- and two-threshold cases 
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Perturbation theory for the Siegert pseudostates (SPS) [Phys. Rev. A58, 2077 (1998) and 
Phys. Rev. A67, 032714 (2003)] is studied for the case of two energetically separated thresholds. 
The perturbation formulas for the one-threshold case are derived as a limiting case whereby we 
reconstruct More's theory for the decaying states [Phys. Rev. A3, 1217 (1971)] and amend an error. 
The perturbation formulas for the two-threshold case have additional terms due to the non-standard 
orthogonality relationship of the Siegert Pseudostates. We apply the theory to a 2-channel model 
problem, and find the rate of convergence of the perturbation expansion should be examined with 
the aide of the variance D = \\E — A™.E'™'[| instead of the real and imaginary parts of the 
perturbation energy individually. 



I. INTRODUCTION 

Resonances occur in a variety of fields of physical sci- 
ences. Despite their diversity, they are characterized 
by two parameters, the resonance energy position and 
width, apart from the coupling with the background con- 
tinuum represented by the Fano profile[l(. A great deal 
of discussions have been given to the interpretation of 
resonance phenomena^]. The most familiar parameter- 
ization of the resonances is condensed into the disper- 
sion formula due to Breit and Wigner. Back in 1939, 
Siegert @] developed a compact mathematical viewpoint 
for characterizing resonances as singular points of the 
dispersion relation. His idea requires the solution of 
the Schrodinger equation subject to the outgoing wave 
boundary condition, 



dr 



ik I <p(r) 



0. 



where a is the radius beyond which the potential en- 
ergy is negligible. The solution 4>{r) is called the Siegert 
state (SS) and it behaves like e zkr near r = a and be- 
yond. This boundary condition destroys the hermitic- 
ity of the Hamiltonian, thus entailing complex-valued 
eigenenergies, i.e., 



E 



2 



E rc 



This is a most direct representation of both the reso- 
nance position and width. This mathematically appeal- 
ing representation had been implemented with tedious 
iterations due to lack of suitable computational tech- 
niques until Tolstikhin et al 3] made a breakthrough 
by introducing Siegert pseudostates (SPS) for the one- 
threshold case. Their idea incorporates the boundary 
condition into the Schrodinger equation so that the dis- 
persion relation is obtained by a single diagonalization 
of the Hamiltonian matrix. Previous applications of SPS 
to resonances in three-body Coulomb problems indicate 
that it is not only a valid procedure but also a new per- 
spective for the SPS representation of resonances and de- 
cay processes Q . Another immediate application of the 



SPS theory is to the time-dependent problem|^,|^] where 
the reflection off the exterior boundary incurs numerical 
instability. Tanabe and Watanabe succeeded in de- 
scribing the reflectionless time propagation based on the 
Siegert pseudostates. Indeed, applied to the half-cycle 
optical pulses, the Siegert boundary condition indeed was 
seen to eliminate the outgoing wave component perfectly. 

Recently, Sitnikov and Tolstikhin ||| Oj stretched the 
scope of the SPS theory by enabling the treatment of the 
two-threshold problem. Despite such progress, there re- 
mains in the theory of SPS a chapter still incompletely 
worked out. This is the Siegert perturbation theory. A 
pioneering work on this subject is due to MorejTfl HU 
who extended the Siegert state theory specifically to han- 
dle the decaying state. The main purpose of this pa- 
per is to complete the Siegert perturbation theory from 
the recently developed SPS viewpoint for both one- and 
two-threshold cases. Particularly, in the one-threshold 
case, we are able to reconstruct More's theory for decay- 
ing states in terms of SPS but with an unexpected 
amendment to his theory. The SPS perturbation the- 
ory (SPSPT) is by no means straightforward owing to the 
non-standard orthonormality of the eigenfunctions. This 
constraint also serves to fix the phase of the perturbed 
wavefunction, a feature which is absent from the stan- 
dard perturbation theory. It is hoped this paper serves 
to expose such noteworthy features of the SPSPT. 

This paper is thus constructed as follows. In SectionlTTl 
we review some basic ideas about the SPS as needed for 
an elementary presentation of the perturbation theory. 
Section ITTT1 gives the details of the SPSPT for both one- 
and two-threshold cases. And Section IIIII deals with a 
specific mathematical model as an example of the SP- 
SPT. Atomic units are used throughout. 



II. THE SIEGERT PSEUDOSTATES 

Since the two-threshold SPS theory contains the one- 
threshold case in itself, we review the two-threshold case 
only, leaving the one-threshold case as the limit where 
the two-thresholds become degenerate 
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A. Mathematical Settings 

Suppose first that there are as many as q independent 
channels. The Schrodinger equation reads 



Id 2 

-2^ + V ^~ E 



where 



V(r) 



<f>(r) = 



v 12 v 2 



\v lq v 2q 



(r) = 0, 



(i) 



V-. 



2q 



Va 



and Vi pertains to the potential energy of channel i, and 
Vij represents the interchannel coupling between chan- 
nels i and j. We consider the situation where there are 
only two energetically distinct thresholds so that we sep- 
arate Vi into two groups. A first group contains channels 
1, . . . , qi and they converge to v\ as r — > a while the other 
group contains channels q± + 1, . . . , q and they converge 
to t>2, that is 



91+1, 



lim V(r) — diag[«i, 



,Vl,V 2 , 



,v 2 ], 



where v\ and v 2 are the two constants representing the 
threshold energies. This allows us to use the 2-channel 
SPS scheme even in the presence of more than two chan- 
nels. The two channel momenta are k\ = \/2(E — vi) 
and k 2 = y/2(E — v 2 ). The boundary conditions are thus 



at r = and 



&(0)=0 



= 



at r = a where j ' = 1 for the first group, i = 1, . . . , q\, 
and j — 2 for the second group, i = q± + 1, . . . , q. Now, 
consider to expand the wavefunction fa by a complete 
orthonormal basis set {717 (r), (I — 1, N)} over r £ [0, a] 
such that 

N 



{ (r) = ^c. M 7r ; (Y). 



1=1 



Substituting this into Eq. Q), and integrating over the 
interval [0, a], we obtain the M = q x A-dimensional 
eigen value problem, 



H - -B - EI M 
2 



c= 0, 



(2) 



where 



H 



( [/(12) 



[7(29) 1 



y u {q2) ■ ■ ■ J 



B = 



1 • • ■ qi qi + 1 • • ■ g— th block 
(k x L \ 



k 2 L 







V 



/ cm \ 



Cl.JV 



Cg,l 



and 



H 



(n) 



1 /" a C?7Tj d7Tj 



dr dr 



■dr 



TTiVnTTjdr 



U 



(mn) 



ITiVmnlTjdr 



■Ki{a)iTj(a). 



In Eq. pj, Zjvf is an M-dimensional unit matrix. The 
eigen system Eq. J2J involves a pair of eigenvalues, A:i 
and ^2 , which may be rewritten as a standard eigenvalue 
equation for a single variable u according to the following 
heuristic procedure. Let us note that energy E can be 
represented by both k\ and k%, namely 



E = -k{ + v 1 



v 2 . 



so that 



where 



{ki + k 2 )(h -k 2 ) -4A 2 



(3) 



A = 



v 2 - Vl 



(Here, we assume v 2 > v\ for simplicity.) Since the prod- 
uct of linearly independent combinations of k\ and k 2 
becomes constant, we require k\ ± k 2 to satisfy the fol- 
lowing conditions, 

k\ + k 2 = 2iAu 
ki-k 2 = -2iAu _1 . 
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Thus, 



and 



with 



fci = iA (u — u 1 ) 
fc 2 = iAiu + u^ 1 ) 



E = v — A 



2 l + » 4 
2u 2 



This procedure of replacing a pair of variables k\ and k 2 
by a single variable u is called uniformization. 



B. The Tolstikhin-Siegert equation 

The uniformization described above reduces Eq. (J2J) to 
M{u)c = Q (4) 

with 

M{u) = I M +uB- +u 2 A + u 3 B + +u 4 I M , (5) 
where 

// : r J1 --/ w r 11 - 1 

A = 



2 

A 2 " 



^(12) H^-vI M ■■■ 



and 



1 



<7i f/i + 1 • • ■ th block 



1 

A 



L 



\ 



By introducing a new vector 



/ 1 

uc 

\ U 3 C 



the non-linear eigenvalue problem, Eq. is reduced to 
a linear one such that 
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u lii 
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2 -> 
U C 
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-Im 


-B- 


-A 


-B+J 




u a c / 




\ W 3 C 



(6) 



Furthermore, the above equation is symmetrizable as fol- 
lows, 



/ 

I M 

I M B- 

\hi B~ A 

( Im 
Im B- 
I m B- A 

V o o o - 



Im n 




/ c \ 
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Im I 




V M 3 C / 



(7) 



Let us refer to Eqs. (QJ, ©, and Q as the Tolstikhin- 
Siegert equations (TSEs). 



III. FIRST AND SECOND ORDER 
PERTURBATION THEORY 

A. Derivation of Perturbation Formulas 

Let us formulate the perturbation theory as appropri- 
ate for the SPS whose orthonormality relation is different 
from the standard one. Relegating the one-threshold case 
to the next subsection, we treat the general two-threshold 
case. We assume the perturbing potential energy V (r) 
vanishes beyond r = a, i.e., 



V'(r) 



(V{ x V{ 2 

v 12 v 22 



KJ 



(r > a). 



The TSE for the n-th state including perturbing potential 
energy reads 



hi +u n B~ 
where 



A 2 



+ u A n B+ + uH M c„ = 0,(8) 



U 



u' 



f(mn) 

ij 



( U' {11) U' {12) 

[/'(12) £/'(22) 



i [/'(I?) [/'(2g) . . . u'{qq) J 



Differentiating Eq. JSJ with respect to A and using the 
orthonormal relationship (see Eq.(44) in Ref. @), 



Im 



u m u„(B - u m u n B + ) 



(u m + u„)(l -u 2 m u 2 m ) 



= Sn 



we obtain the Hellmann-Feynman theorem (HFT) in the 
present context, namely, 



c*U'c n = A : 



1 



,1 + ut 



= u_ A 2 : 

dX dX \ 2u 2 



dE n 
dX 



(9) 
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Now, we consider the perturbation series of u n and c„ 
such that 

u n = v,™ + \u£ } + \ 2 u^ + . . . , (10) 
c„ = c <o) + A c < n 1 )+A 2 cl 2) + ---, (11) 

where uffl and ci "* are the n-th solution to the unper- 
turbed equation, Eq. J2J, 

M(4°))cf ) = 0. 

Substituting the perturbation series, Eqs. (fTU|) and lfTT^ I. 
into Eq. and then comparing each power of A, we 
obtain 



A : A ; 



A 1 



1 — u 



(0)4 



,(0)3 



(12) 



(0)3 



2« 



2^(1-^)- 



(13) 



Next, let us evaluate the expansion coefficients over the 
unperturbed eigenstates. To this end, we rewrite Eq. © 
using Eq. (0, namely, 



M(u n )cn 



U C " 



so that 



2\ulM-\u n ) T . 



U'c n . 



The spectral representation of Ai is given by 



4M 



M-\u n )=Y, 



a i c i c i 



^2{l-uf^){uP-u n ) 
(See Eq. (59) in Ref. 8].) Using the relations 



4M „(0W0)-<0)T 



1=1 1 - M ( 



(see Eqs. (51) and (52) in Ref. @|), we have 



iM 



\M- 1 (u n )=Y, 



(0)3J0)J0)T 



(14) 



(15) 



(16) 



Substituting this into Eq. (|14f> and comparing both hand 
sides power by power for A, and then using Eqs. I|12|) and 
(fHIjl . we have 



A 1 : & 



1 „(0)3aC0)-s(0)T 
4M (0)3 



c<°) 



A2 E ,. |") („!,'/ 

A (I""/ ){v> n > -U\ ') 

A2 X/ 



d 0) 



,(0)3 



Aa fc(i-«r)(«g ,) -«[°>) 1 



(17) 



,(0)4 



24 0) (i-ur" i ) 



i 7/ (i)^o) 

(0)4s " " 



-cf > + 



nn 



l 



l 



(0)2 ,(0)2 



1.(0)1.(0) 



do) 



(18) 



where 



W" = ^ 0)T C/'c< 0) 

r r mn rn 



and, as before, 
fc^^zA^)-^))- 1 ], 



Let us note that for , there is a term on top of the 
summation, which is made absent in a standard pertur- 



bation theory because the normalization is unchanged in 
so far as this term is purely imaginary under the standard 
orthogonality relation. This freedom is not warranted in 
the present case. 



Finally, we have the perturbation formulas for the two- 
threshold SPS, 



J 



Eii) = 40)^0) 



(19) 



5 



(2) 



AM 

W' 2 ( 1 1 



v m W' 2 

U l VV ln 



(1 



,(°) 4 v„,(°) 



',(0)2 
'l/i 



u(0)2 
v 2» 



1.(0)1.(0) 



(20) 



B. One-threshold case as a degenerate limit 

It is important to clarify the relationship between one- 
and two-threshold cases. In the following, we prove that 
perturbation formulas for the one-threshold case are ob- 
tained when we implement a limit of vi — > v\. In this 
limit, the following scaling clarified in Ref. @, 



reduces the two-threshold TSE to a one-threshold one, 
namely, 



(A + kB + k 2 I m )c = 



(22) 



c 



a 





where 



ra.u 



(12) 



.4 



/ #"« - vI M U^ 1 

u{ i2) hW-vIm 



t/(2g) 



f/(2«) 1 



• • • #W - vI M J 



FIG. 1: Broken lines: Diagonal elements of the potential ma- 
trix in Eq. Solid lines: adiabatic potential energies. This 
system supports three resonances: shape type (a) in channel 
1, Feshbach type (b) and shape type (c) in channel 2. 



and 



B 



1 • • • q— th block 
-L ■■■ 

••• —L 



and k = ik\ = ik^. This scaling corresponds to the so- 
lution k\ = &2 in Eq. © when v\ — > i>2- Note that the 
solution fci = — hi in Eq. @ is unphysical since the de- 
generate threshold here means the equivalence of asymp- 
totic wavefunctions in this limit. 

Thus, the scaling leads us to the perturbation formulas 
for the one-threshold case, namely 



2M 



E {i) 

n 



W 

nn 

2M 



k (0) (k {0) 

ft,] \rvn 



So) 



-Ci 



W nn -<0) 
9i .(0)2 C « 



w{ 2 



Ijtn ft Z 



(fcr-fc ; (0) ) 



W' 2 

r r nn 

oju(O) 2 ' 



formulas |10|. Our expressions for the first-order eigen- 
vector and for the second-order eigenenergy are differ- 
ent from his|ll|. The origin of the discrepancy has been 
traced to an algebraic error in More's derivation of the 
first-order wavefunction. (One necessary term is unfor- 
tunately dropped during his derivation.) As a result of 
this, an extra term is restored in either formula. Here, 
one important difference from the standard perturbation 
theory is that no Hermitian conjugates appear in these 
formulas. This might suggest at first that there would 
remain phase ambiguity. However, any ad hoc additive 
phase would instead mar the orthogonality relation, that 
is what is the relative phase in the standard theory is 
fixed in the SPS theory, thus leaving no ambiguity with 
the phase of eigen functions. It is thus worthwhile to 
see the consistency of the orthonormality relation and 
the Siegert boundary condition for the particular case of 



#) 



This verification is worked out in Appendix. 



A Model Problem 



Note that the summation runs over the branch of k\ = ki , 
that is only over a half of the full non-degenerate space. 
These correspond to the SPS representation of More's 



Let us present an example of the perturbation the- 
ory for the two-threshold case. We revisit the 2-channel 
model potential with two thresholds that is taken up in 



6 



Ref. 



i.e. 



V(r) 



15e -0.5r 

5re _r 



5re r 
15(r 2 — r — l)e~ 



15 



(23) 



The potential V(r) supports three resonances. The adia- 
batic potential energy curve of the first channel supports 
one shape type resonance (a) while the other channel 
supports one Feshbach type (b) and one shape type (c) 
resonance. These resonances are depicted in Fig. ^ Wc 
carried out the diagonalization of the TSE, Eq. (JBJ, using 
the discrete variable representation (DVR) functions as a 
basis set. The calculated resonance energies and widths 
with different numbers of the basis functions are given 
in Table [I] Let us call these results as direct numerical 
solutions. To implement perturbation calculations, we 
separate V{r) into 



V(r) = V Q (r) + V (r) 



where 
Vo(r) 
V'(r) 



15e -0.5r 4re -r 

4re' r 15(r 2 - r - l)e~ r + 15 

re~ 
re~ r 



perturbation energies using the unperturbed solutions of 
TSE for the same box size a = 50 as in Ref. 8j. Table 
U shows the results of first- and second-order perturba- 
tion calculations, and Figs. 12141 depict how the numerical 
solutions converge in the complex plane. In the present 
model problem, the first-order resonance energy agrees 
with the direct numerical solutions to about 2 to 4 dig- 
its while the width agrees to about 2 to 3 digits. And 
the second-oder resonance energy agrees to about 3 to 
5 digits while the width agrees to about 1 to 3 digits. 
An important fact which we must remark is that the res- 
onance energy and width do not appear to converge in 
pace. For instance, the width of resonance "c" evaluated 
by the second-order perturbation theory appears less ac- 
curate than the first-order one while the resonance energy 
appears to have improved. The seeming deterioration of 
the width is a little overwhelming, all the more so for the 
improvement of the resonance energy. Nonetheless, the 
distance between the second-order result and the direct 
numerical one becomes rather small (see Fig. 0J) in the 
complex plane, that is in the Siegert state perturbation 
theory the convergence is to be measured with respect to 
the variance 



D = \\E-^2x n E^ 



(24) 



We regard Vq as the unperturbed potential energy and 
V as the perturbation potential energy. We calculate 



rather than with respect to the real and imaginary parts 
of the sum, individually. 



IV. CONCLUSIONS 
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APPENDIX A: CONSISTENCY WITH 
ORTHONORMALITY RELATIONSHIP AND 
SIEGERT BOUNDARY CONDITION IN FIRST 
ORDER 



It is of interest to speculate on possible uses of SPSPT. 
One immediate application would be to the manipulation 
of Siegert poles. The shadow poles located near the phys- 
ical sheet may be transformed to physical resonances by 
an appropriate perturbation. We leave issues such as this 
for a future task. 



Here, we prove that the first-order wavefunction sat- 
isfies the orthonormality relationship and the Siegert 
boundary condition consistently. First of all, we expand 



■N 



-B 



(Al) 



7 



TABLE I: Columns 5R,. 9, D, and ^represent the real and imaginary parts of resonance energies, error variance in the complex 
plane, and the dimension of the DVR basis set, respectively. 
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_E (Direct numerical solution) 
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FIG. 3: Complex energies for resonance b 
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